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1 Introduction. 

Markovian random evolutions in R n were studied in the work [1], where the connec- 
tion between the equation for the functional of the evolution and Bessel equation of 
high order was found. 

Bessel equations of high order and their solutions - Bessel functions of high order 
- were studied by Turbin and Plotkin in [2]. 

In the works [3,4] Orsingher and Sommella managed to combine these two results 
and to receive the obvious view for the distribution functions of random evolutions 
in R 2 and i? 3 . In [4] they also made a conjuncture about the view of the distribution 
function of the random evolution in R n . 

In this paper we generalize the results of Orsingher and Sommella and prove their 
conjuncture. 

Let us start with the definition of Markovian random evolution. To do this we 
have to define a regular n + 1-hedron inscribed into a unit sphere in R n . 

Definition 1.1: We call figure in R n a regular n + 1-hedron inscribed into a unit 
sphere if the following conditions are true: 
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1. The figure has n + 1 vertices tq, . . . ,r n situated at the unit sphere with the 
center in (0, . . . , 0). 

2. The center of masses of the figure is situated in (0, . . . , 0). 

3. The i-th vertice is situated into i + 1-dimensional subspace, determined by the 
first i + 1 coordinates of R n . 

Having this definition, we may easily find the coordinates of the vertices. The 
components of Tj, i = 0, n are equal to 



Really, the first condition is true, because the distance from every vertice to the 
zero point is equal to 1. 

The second condition is true, because the sum of first components of all coordi- 
nates is equal to 0, the same for second components, etc. 

At last, for the i-th vertice i + 2, i + 3, . . . components of coordinates are equal to 
0. So the third condition is also true. 

We are now ready to define Markovian random evolution in R n from [1]. 

Definition 1.2: We call a random process S(t) Markovian random evolution in 



1. The process begins its motion at the point x = (xi, . . . , x n ). 

2. The initial direction of the motion is Tj,i = 0, n, where the components of 
coordinates of Tj are presented in (1). 

3. The time of the motion at some direction is distributed like e~ xt . 

4. The A;-th direction is followed by the k + 1-th direction (the ro-th one is followed 



5. The velocity of the particle's motion is equal to v. 

The system of backward Kolmogorov equations and corresponding n+l-dimentional 
equation for the functionals (here % = 0, n is the start direction) 



were also studied in [1]. 

Using a well-known result from the theory of Markovian processes we may receive 
the system of differential equations for the distribution functions 

fj(xi, . . .,x n ,t)dxi . . . dx n = P{Si(t) £ dxi, . . . , S n (t) £ dx n ,D(t) = j}, (3) 

where D(t) is the direction at time t. 

The system of equations for the distribution functions (3) is adjoint to the system 
for the functionals (2) (see [1]): 




(1) 



R n if: 



byO). 



Ui (x,t) =Ei(f(S(t)) 



(2) 



'&/oCM) 



A/o(M) 
Xh(x,t) 



v{T0, V)/o(x, t) + Xf n (x, t) 

v{T 1 ,V)f 1 (x,t) + Xf (x,t) 



(4) 




Xfn(x,t) 



v(T n ,V)f n (x,t) + A/ n _i(z,f) 
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The corresponding n + 1-dimentional equation is: 

II {j t + A + v{r h V)) /(x, t) = X n+1 f(x, t), 



(5) 



where / is any of functions fo, ■ ■ ■ , f n - 

To describe the support of the distribution of random evolution we use the notion 
of regular n + 1-hedron from Definition 1.1. 

One of the most important characteristics of random evolution that differs it 
from Brownian motion is that evolution with probability 1 moves into a regular 
n + 1-hedron 



<-vt 



vt 1 

xi,...,x n : < xi < vt, — — 

n n — k + 1 



(n-1)... (n-k + 1) 
(n + 1) . . . (n - k + 3) 



x±+ 



(n-k + 1) (n-l)...(n-k + l) 



(n-A; + 3)' 



(n + 1) . . . (n - fc + 3) 



< 



(n-l)...(n-fc + l) (n-k + 1) 



(n + l)...(n-fc + 3)" 



(n-fc + 3)' 



^/(■.-l)-(n-*+l) iit = ^ 



(n + l)...(n-fc + 3) 
The probability that evolution is on the edge dT vt of n + 1-hedron is equal to: 



(6) 



P{S(t) £ dT vt } = e~ xt + Xte~ M + ...+ 



-At 



(Xt) 



n—1 



-Xt 



(n-1)! 



Here e is the probability of being at time t on the vertices of T vt , Q^jj) ° 
1, n — 1 is the probability of being on /c-dimensional edge of T vt . 

So, the continuous part of the distribution we are to find has the property: 

(xt) n - 1 



J ... J /(xi, . . .,x n ;i)dxi . . . dx n = 1 - e 



xt - Xte~ xt 



(n-1)! 



(7) 



In other words, we have to find a non-negative continuous function that satisfies 
system (4) and condition (7). 



2 Derivation of the equations. 

Making exponential transformation fj = e~ xt gj,j = 0, n we receive from (4): 
M9o(x,t) = -v(T ,V)g (x,t) + Xg n (x,t) 



d - t gi(x,t) = -v(T 1 ,\7)g 1 (x,t) + Xg (x,t) 



(8) 



m9n(x,t) = -v(T n ,\7)g n (x,t) + Xg n -i(x,t) 



\ dt 



Theorem 2.1 The functions gj,j = 0, n are solutions of the equation 

n+l / fz: \ n+l / i \ n+l 



dg 



dyi--. dy n+1 



>n 



where 



Vk 



Vn+l 



vt 

yi = — \-xi, 

n 



1 



9, 



(9) 



n-k + 1 



fc + 1) 

(n + 1) . . . (n - k + 3) 



(n-l)...(n-fc + l) 
(n + 1) . . . (n - + 3) 



n-k + 1 
n — k + 3 



+ x k , k = 2,n, 



Vn+l 



'(n-l)...l 
(n+l)... 3 



'(n-l)...l 
(n + l)... 3 



xi - . . . - 



2^n— 1 X n . 



(10) 



Proof. By applying the transformation (10) to the system (8) we receive: 

n+l dgp _ \ 



v 7 " 



yn+i 



1- 



By differentiation we easily have (9). 
Theorem is proved. 

Theorem 2.2 The transformation z = n+ ^/yi ■ ■ ■ y n +iconverts the equation (9) 
into n + 1-dimensional Bessel equation: 



d_ 

dz 



n+l 



9 



n+l 



(y/n{n+ 1)) 



■\ n+l 



1 



n+l 



(11) 



where g is any of the functions gj,j = 0, n. 

Proof. We don't show the proof since plain calculations are involved. 

The solutions of equation (11) were found in [2]. We use here one of the functions 
that is solution of (11): 



'0,n 



A Vn(n + 1) _ 
v 2n+ %2n + 2 



E 

fc=0 



A v / n(n+ 1) 1 
v 2 ™V2n+^n + l' 



(n+l)fe 



(fc!)™+ 1 ' 



(12) 



Remark 2.1 The proofs that this function is solution of (11) in case of n = 2,3 
may also be found in [3,4]. 

Making the backward change of variables in (12) and exponential transformation 
we receive the solution of system (4). 
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3 Probability distribution. 



Now we are ready to formulate the main result of the article: 

Theorem 3.1 The absolutely continuous component of the distribution of Marko- 
vian random evolution in R n is: 



f(xi, • • • , x n ] t) 



n„— At 



(y/n+~T) n+1 v n 



\ n + \ r - 



-1± 

at 



+ ...+ 



+ 



Qn 



A ^n(n + 1) 
v 2n +^2n + 2 



^Vi ■ ■ ■ Vn+l , 



(13) 



here and later i = 1, n + 1 are the linear combinations of t, Xj, z = 1, n pointed in 
(10). 

Proof. The fact that (13) is a non-negative function may be easily seen from the 
view of the series (12). 

The function (13) satisfies the equation (5). Really, 7 ,n (j; ^+^^ x x n+ \/Vi ■ ■ ■ Vn+i) , 

where yi,i = l,n + 1 are pointed in (10) satisfies equation (9). Any linear combina- 



tion of derivatives of the last function with the coefficient 



also satisfies 



v /n+l) n + 1 v n 

(9). So, if we make backward exponential transformation, we receive a solution of 
(5). 

The last problem is to prove the condition (7) for the function /(xi, . . . , x n ; t). 
We first note that 

f f T f\ Jn(n + 1) , \ , 

J -J hn 1^-1^ "WTT^T j dx 1 . . . dx n = 



E 

fc=0 



^ N (n+l)k 



(n + 1)^7 ( 2n+ v / 2n + 2)(™+ 1 ) fe (A; 

[ V ["+^:::4 a: i+-+v / l a: "- 2 ~ v 7 [^H H 



1 [ / (n-l)...: 

2 LV (n + l).., 



/ 



1 

2^ 



-\l -X n - 2 + 



/(n-l). -.2 



'(n- 1)...2 
(n + 1)...4 



Xl 



(n + l)... 4 

/ 

f /(n-l)...l /T / (n-l)...l ,1 

[V U+l)-3 3:i + - + V S^-^V (n+l)-3 ^J 



ttf + 2x n _i 



dx n _ix 



'(n-l)...l 
(n + l)... 3 



xi- 
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- \l-x n -i + 



(n + l)...3 



— x^ 



dx r , . 



(14) 



The inner, first, integral could be found using the change of variables x n 

(-VferH Xl " • • • " ^K-i + Vfew3 vt ) z - We have: 

/ (n-l)...l . /T / (n-l)...l ,1 

V (u+l)...3 a:i + - + V 3 X "-1-V U+1)-3 ^J 



+ 1J...3 



Xl+. 



3 1 "- 1 V (n+1). ..3™ 



! {n- l)...l 

(n + 1)... 3 



Xl 



'(n-l)...l 
(n + l)...3 



x r . 



dx r 



'(n- 

(n + 1). ..3 



Xl 



3 X "- r 



Zn-i+ 



! (n- 

(n + l)...3 



)2fe+l 1 , 
y (i-^) fc dz= f- 



+ W -x„_i + 



(n-l)...l 
(n+1)... 3 



'(n-l)...l 
(n + 1)... 3 

> 2fc+l 



>2jfc+l 



ft 



(*!)' 



(2fc + l)! 



In the following, second, integral 



I 



i 

2fe 



l[ / (n-1)... 
2 [V (n+1)... 



Il + ...+ i/|l„-2' 



/> 

"V 



+ 1J...4 



tit 



'(n-1)... 2 
(n + 1)... 4 



Xl 



+ 1 



+ 1 



'(n-1) 



(n + 1) 



-ui + 2x n _i 



'(n-l)...l 
(n + 1)... 3 



2fe+l 



r [WW- 



dx n -i = 



(n-1)... I 
(n+1)... 3 

V 

/ 



Xl 



-x n -\+ 



-^i+-+A/-^-2- A / ; n+1 ;... z 



i )- 2 3;i i + ./i x 



2 fc \/3' 



■2fe+l 



'(n- 1)...2 
(n + 1)... 4' 



xi - . . . - 



+2x n _i]' 



•Era— l] dx n ^i 



(n-1)... 2 
(n + 1)...4 

1 



Xl 



Zn-2 + 



'(n-l)...2 
(n + 1)... 4 



vt+ 



'(n-1) 



'(n-1)... 2 
(n + 1)... 4 

^2 



2 fc \/3' 



2fc+l 



(n+1) 



-xi - . . . - 



X n -2 + 



+1 



, v 3fe+i - y t+ 1 i)':.A xi+ - + ^ Xn - 2 - V 



(n + l)...4 / 



/ 



(1+ 



+2 



%n—l 



■l)-2 



+l)..-4 



).</. 



(n-l)-2 
(n+l)...4 X l 



2fc+l 



x 1- 



33fc+l 



%n—l 



(n-l)...2 
(n+1). ..4 



Zi - . . . 



4 ^n-2 + y ( n+1 ) 



(n-l). ..2 



ft 



'(n-l)...2 



-xi - . . . 



2 3fc+i^/3 2fc+1 y y (n + 1)... 4 

/ 



_ / x 3fc+l 

2 /(n-l)...2 \ 



1 2 
- + -x 
3 3 



%n—l 



(n-l)...2 /2 , /(n-l)...2 , 

Wrpi Xl - ■ ■ • - V 4 X ™" 2 + v teb " f 



2 2 

3 ~ 3 



£n-i 



(n-l). ..2 
(n+1). ..4 



4^-n— 2 



we make the change of variables z - 
so we have 

3 3fc + 2 / /(n-l)...2 



1.2 

3^3 



(n-l)...2 
(n+1). ..4 

x n — 1 




2fc+l 



n-l 



(->/§ 



-l)-2 



+ 1J...4 



II- 



3fc+2 



and 



2 3fc+2^/3 2fc+1 ^ Y (n + 1)... 4 
i 

33fc+2 



12 /(n-l)... 2 \ 

4 X - 2 + V (n + l)...4 ^ 



y z k {i - z ) 2k+i dz = 



(n-l)... 2 



2 3fe+2 v /3 2fc+1 I Y (" + !)• --4 

3fc+2 



Xi- 



/ (n-l) ...2 ^ r(fc + l)r(2fc + 2) 
|(n + l)...rj T(3A: + 3) ' 

The following integrals could be found in the same way. For the m-th integral we 
have: 



m — 1 
m + 1 



m(fc+l)-l 



m + 1 \ ( m + 1 )( k + 1 )' 1 Y(m(k + l))r(fc + 1) 
m J r((m+ l)(fc + 1)) 




(n — 1) . . . m I m 

%1 • • • \ / ; 7T%n—m~\~ 



(n + l)...(m + 2) "' V m + 2 

, v (m+l)(fc+l)-l 

/ (n — 1) . . . m \ 

+ V (n+l)...(m + 2) Vt ) 

Substituting the integrals found into (14) and performing simple calculations we 
receive: 




A -von = ar { -^^x 



x 



~ (At) (n+l)(fc+l)-l 

^((n + l)(fc + 



But to prove (7) we have to find the following integrals: 

r r t) m t I \ \/n(n+l) \ . , q 

J • • • Jr„ t aF^ J 0,n U ^+^ 2ra +2 " + ^ T7 ^+1 J dxi . . . dx n , 171 = 1,1%. 

To do this let us find 

3 /" f T f X Jn(n + 1) ; \ , , 9 f 

diJ -J^Av ^2n + 2 " + ^ TTT ^ )dxi...dx n = -J ... 



I 



h,n ( - X 



\fn{n + 1) 



n+ \/yi • • • y^+i H Xi • • • dx n -idx n . 



The main rule for such integrals is: we find the derivative of the above limit and 
multiply it on the integral in which we substitute the above limit instead of x±, then 
we subtract the derivative of the down limit multiplied on the integral in which we 
substitute the down limit instead of xi.On so on for every integral. 

But in our case when we substitute the above limit into the next integral it 
becomes equal to - this may be easily seen using (6). Otherwise, when we substitute 
the down limit the argument of Bessel function becomes equal to - this may be seen 
from (6) and (10). But it's clear that io,n(0) = 1, so the integrand will be always 
equal to 1. 

In other words, we have the following equality: 



d_ 

at 



f f r f X Jn(n + 1) , \ , 

]-]h,n y- " WTTT^ j dXl ... dx n 
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d f /" , f f d r f X Jn(n + 1) , \ , 

= -J ... Jld Xl ...dx n+ J ... J -I , n I - ^Vi ■ ■ ■ Vn + 1 U Xl ... dx r 

T vt T vt V 7 

But the first integral is the derivative of VolT vt . We may easily find that VolT vt = 
So, 

a /• [ T f x V n ( n + 1) , \ , 



dt 



Un+T) 



For the next derivative we have: 



I f In ( x , +] 

y - y io '™u 2 »v2^+2 v 



y/yi ■ ■ ■ \dxi . . . dx n 



(V^+1) W+1 v n f n-2 | g 

(V^) n (n-2)! dt 



a^ + 2 dXl ■ ■ ■ dXr ' 



We'll find the last integral by the rule we mentioned above, but here the above 
limit makes the following integral equal to 0, and when we substitute the down limit 
into the integrand we have |/o,n(0) = 0. So, 



d 2 f f I X Jnin + 1) ; \ , 

w J ■ ■ ■ J Io > n [v + 2 b+ ^ 1 -^ 1 ) dxi • • • Ax ' 



_ v n t n-2 + 



(y/n) n (n - 2)! 

'J- vt 

By analogy, for the k-th derivative we have (k = 3, n): 

X Jn{n + l) 

I J 0,n I 

(A+T) n+1 



d 2 T (X Jnin + I) , \ , 

dT 2 ^ n v + 2 j dxi ■ ■ ■ dXn - 



T 7 . V 



2 "V2^+2 



n+ y/yi ■ ■ ■ Vn+1 \dx 1 ... dx n 



{y/n) n {n - k)\ 



v n t n-k + 



f f 9 

J ■■■y w 



• k T fX Jn(n+1) , \ , 

J 0,n I ~ n+ V / yi • • • ^1 • • ■ ^r. 



v 2n+ v/2n + 2 



Using the formulas found, we receive: 



5t 



'A Jn(n + 1) 



v 2n+ v /2n + 2 



• i+ v / yi • • • H Xl • ■ ■ dXn = 



9 



= V ' 



-xt 



- 1 - At - ... - 



(At)"' 1 ' 
(^T)!_ 



And finally, 



-At 



(VnTT)"+ 1 t; n 



A™ + A n ~ 



at 



... + 



gn 



Y / n(n+T) 



• • • H Xi • • • dx « = 1 ~ e 



-At 



r ■ ^ 

(At)"' 1 

(n-1)!* 



-At 



Theorem is proved. 

Remark 3.1 In the work [4] Orsingher and Sommella conjectured that the normal- 



-. As we've seen it is not 



izing constant in (13) should be equal to — r - ■ - 

° v ' 1 (vn+l) n+ L v n nWollvt 

surprising because this constant appears when we find the integral f ■ ■ ■ J T ldx\ . . . dx r , 
VolT vt . 
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